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Abstract 


Parasitic plants in Scrophnlariaecac and Orobanchaeeae have been traditionally depicted as forming a linear evo¬ 
lutionary series beginning with herniparasitism and ending with holoparasitism. The genera Ijithraea . Harveya* and 
Ityobanche have been viewed as transitional links between the parasitic members of Scrophulariaceae and the strictly 
holoparasitic habit of the traditional Orobanchaeeae. Phylogenetic analyses of plastid rps 2 and null K sequences were 
performed. Ihe transitional genera are not transitional to the traditional Orobanchaeeae, but represent multiple inde¬ 
pendent origins of holoparasitism. W ithin Scrophulariaceae. the two traditional subfamilies Rhinanthoideae and Antir- 
rhinoideae are defined by the arrangement of the corolla lobes during aestivation. However, neither of the two subfam¬ 
ilies is monophvletic in our analyses, suggesting that corolla lobe position is a homoplastic character. While the 
traditional Orobanchaeeae are monophyletic, tribes Buchnereae and Khinantheae are clearly not, and genus Orobanche 
probably is not. Clades of parasitic genera correspond well with biogeographic provinces. One strongly supported clade 
contains the parasitic Scrophulariaceae, the traditional Orobanchaeeae, and Lindenbergia. It is proposed that this clade 
be defined as the Orobanchaeeae. 


Parasitic angiospemis are found in 16 families 


encountered. The group with the greatest range in 
and live in diverse habitats, ranging from tropical parasitic specialization is found within the Scro- 
forests to arctic islands (Mussehnan & Press, phulariaceae/Orobanchaceae. This family pair in- 


1995). Recent work on several parasitic groups has eludes nonparasites, both facultative and obligate 

explored their anatomy, physiology, ecology, and hemiparasites (photosynthetic plants that obtain 

molecular biology and the control of economically water, nutrients, or photosynthates from their host), 

significant parasites (Press & Graves, 1995). Par- and holoparasites (non-photosynthetic, obtaining all 

asites exhibit dramatic adaptations. Some lack photosynthates from their host). Much has been 

leaves, stems, roots, and the ability to photosynthe- learned about the anatomy and physiology of par- 

size. For many lineages of non-photosynthetic par- asitism in this group, especially in Striga Loureiro 

asites, study of their evolution is impeded by the and Orobanche L., which are important pests of 

very fact that they are so specialized and so differ- crop plants (Press & Graves, 1995). Reliable in- 


ent from their nearest photosynthetic relatives. In formation on phylogenetic relationships among the 
many cases the identity of those relatives 



various nonparasites, hemiparasites, and holopar- 
unknown (Nickrent & Duff, 1996; Nickrent et al., asites would enable us to better understand the an- 
1998). The evolution of parasitism can best be atomical, physiological, and genetic changes that 
studied in a group where near relatives are known, occur during the evolution of parasitism and the 
and different degrees of parasitic specialization are loss of photosynthetic ability. 
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Boeshore (1920: 171) argued for a close rela¬ 
tionship Between the Serophulariaceae (figworts) 
and the traditional Orohanchaceae (broomrapes), 
describing in detail an evolutionary transition se- 


“All macroscopic* and microscopic details suggest that 
the parasitic Serophulariaceae and Orohanchaceae form 
a continuous and parasitically degrading morphological 
series that show transitional steps from green nearly au¬ 
totrophic plants like Melampyrum, Rhinanthus , and Eu¬ 
phrasia to increasingly condensed and degraded genera 
like Bartsia and Harveya , on to Lathraea , that has been 

m/ 

shown to be placed by some botanists in Serophulari- 
aeeae, by others in Orohanchaceae, thence through spe¬ 
cies ol Orobanche to Epiphegus , and finally Aphyllon 
and Conopholis . M 


According to Boeshore, the evolutionary series 
among parasitic figworts and broomrapes involved 
progressive enlargement and consolidation of haus- 
toria (the connection to the host), shortening ol the 
vegetative stem, reduction of leaves, simplification 
of the inflorescence, change from few, large seeds 
to numerous, small seeds, and a reduction in the 
ovary from two locules to one. Weber constructed a 

j 

different evolutionary series based on haustorial 
anatomy (Weber, 1980), but otherwise, Boeshore’s 
depiction of the group’s evolution has endured re¬ 
markably well to the present. Later authors have 
cited Boeshore and reiterated his ideas. The para¬ 
sitic Serophulariaceae were viewed as '"a pointer to 
. . . Orohanchaceae'* (Hutchinson, 1969: 630); 
Kuijt (1969: 95) observed, “Many features of Oro- 

banchaceae appear to be already foreshadowed in 
their figwort relatives," and Cronquist (1981: 940) 
stated, “The evolutionary journey toward parasitism 
obviously begins in Serophulariaceae; the Oroban- 
chaceae merely occupy the house at the end of the 


road." 


Recent results from DNA sequencing show a 


more complicated history for the group. Phyloge- 
ies based on the plastid genes rps2 and rbcL 


(dePamphilis et al., 1997; Wolfe & dePamphilis, 
1998; Nickrent et al., 1998) show that parasitism 
evolved a single time in the figwort/broomrape 


family pair. Thus all the hemi- and holoparasites, 
taken together, form a monophyletic group. This 
agrees with Boeshore, whose evolutionary j 


also contains a single origin of parasitism. Within 


this parasitic clade, however, evolution has not 
proceeded as Boeshore envisioned. The rps2 phy- 
logeny also shows that the transitional genera 
( Lathraea , Harveya , Hyobanche) are not transition¬ 
al at all, but occupy their own branches among 
the parasitic figworts, representing separate losses 
of photosynthesis, independent of the broomrap 


(dePamphilis et al., 1997). It is not surprising that 


Boeshore had developed a different view of the 
groups evolution based on morphology. Many of 
the morphological characters that bear on phylog- 
eny are probably adaptations to the parasitic life¬ 
style, and most involve reduction or loss of struc¬ 
tures. This makes it difficult to interpret 
relationships within the group. It seems that, as 
the separate lineages independently specialized as 
parasites, they converged (sensu Eldredge & Cra- 
craft, 1980) on similar morphologies. This paper 
includes additional DNA sequence data, enabling 
us to more rigorously test whether the parasites 
form a monophyletic group and whether the hol¬ 
oparasites form multiple groups. Because the five 
basal branches of the parasitic clade were left un¬ 
resolved in the rps2 study (dePamphilis et al., 
1997), we have, in this study, included a wider 
sample of taxa lor rps2 and added characters from 
another plastid gene, matK. This gives us the in¬ 
crease in resolution needed to provide a more de¬ 
tailed phylogeny. 

A new definition of Orohanchaceae will be pre¬ 
sented. For that reason, the common name 
“broomrapes” will be used for the traditional Oro- 
banchaceae. 


Throughout this paper we w in use the class ifi- 
cation system of Wettstein (1897), which, although 
dated (Thieret, 1967), is the best worldwide treat¬ 


ment of Serophulariaceae. There are two minor ex¬ 
ceptions to our use of Wettstein’s system. We have 
included his subfamily Pseudosolaneae in the An- 
tirrhinoideae (Thieret, 1967), and we use 1 the name 
Buchnereae Benth. for his tribe Gerardieae Benth., 

following Pennells (1935: 379) “desire that the 

tribal designation be based on a generic name un¬ 
questionably applying to a member of this group." 
Serophulariaceae are divided into two subfamilies. 


Antirrhinoideae and Rhinanthoideae, “which differ 
fundamentally in the single character of the former 
having the posterior (corolla) lobes, and the latter 
the anterior lobes external in aestivation." (Pennell, 
1935: 40). This corolla aestivation character has 
been studied throughout floral development for a 
number of figwort genera (Armstrong & Douglas, 
1989), and the ontogenetic* patterns were consistent 
within each subfamily. However, a study of rbc L 
and ndh F sequences (Olmstead & Reeves, 1995) 
found the Antirrhinoideae to be polyphyletic. In 
that study, the Rhinanthoideae, which consist of 
two nonparasitic* tribes (Digitaleae and Veroniceae) 


1 two parasitic tribes [Buchnereae and Rhinan- 
theae Benth. (= Euphrasiae Benth.)], were repre¬ 
sented only by two nonparasites. Digitalis and Ve¬ 
ronica. Thus, neither the position of the parasites 
nor the monophyly of the Rhinanthoideae could be 
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addressed. Phylogenies of r/ts'2 (dePamphilis et al., suit*, and Iwo hractlets subtending the flower (Pen- 
1997, and unpublished) and r/wL (Wolfe & de- nell, 1935). 


Painphilis, 1998) indicate that the parasitic and These questions were addressed using two plas- 
nonparasitic Rhinanthoideae each form a clade, tid genes: rps2 and mat K. These are two of tin 


and the two clades are separated by many other small set of plastid genes that are found intact in 


taxa including Bignoniaceae and Lamiaceae. Co- all ol the taxa in this group, even in the holopar- 
rolla aestivation pattern is also used to separate the asites. Although the rihosomal protein gene rps2 is 
two parasitic tribes of subfamily Rhinanthoideae. new to plant phylogenetic analysis (dePamphilis et 
Of the three ahaxial lobes, the center one folds over al., 1997), matK has become widely used and has 
the two outer ones in the Buchnereae, and the outer many advantages lor use in phylogenetic studies 


two fold over the center one in the Rhinantheae (Hilu & Liang, 1997). It is large (~ 1500 bp) and 

(Armstrong & Douglas, 1989). The membership of rapidly evolving (Wolfe et al., 1992; Olmstead & 

tilt *se tribes has been stable since their inception. Palmer, 1994), changing 2 or 3 times as fast as r/>cL 
The broomrapes or traditional Orobanchaceae (Steele & Vilgalys, 1994; Johnson & Soltis, 1995), 


(Beck von Mannagetta, 1930) consist entirely of thus providing many substitution characters. The 
holoparasites and have long been hypothesized to changes in mtifK are distributed more evenly 
have affinities either with Serophulariaceae or the among the three codon positions and throughout its 


Ges 


eae. Gesneriaceae have been considered length than in other, more conserved genes (John- 

the closest relatives of tin* broomrapes by some au- son & Soltis, 1994; Hilu & Liang. 1997). As a cod- 


thors (Baillon, 1891; Fritsch, 1895; Wettstein, ing region, matK. is 


to align than non-coding 


1897; Beck von Mannagetta, 1930), emphasizing regions, yet most data sets have a few small indels 
the unilocular placentation of both families. Hallier (insertions or deletions) of 3 to 15 bp, providing 
(1903) and Bellini (1907) both proposed formal in additional characters. Int erestingly, although the 


elusion of the broomrapes in subfamily Rhinan- 




thoideae of Serophulariaceae. Boeshore (1920) also 1991), and possibly several other Orobanchace; 


found the broomrapes to be similar to the parasitic taxa (N. D. Young & C. W. dePamphilis, unpub- 
figworts, not only in the parasitic habit and the lished), we found an intact matK gene in all Oro- 


presence ol haustoria, but in many other features, banchaceae and Serophulariaceae tested, including 


such as: (1) plants similarly reduced, with very 



sue 


h 




(dePamphilis & 


short vegetative axes and linear inflorescence axes; Palmer, 1990; Wolfe et al., 1992) and ConopholLs 
(2) leaves reduced to tooth-like scales; (3) sepals (Colwell, 1994). 
and stamens with tapered, multicellular hairs in¬ 
terspersed with capitate-glandular hairs; (4) paral- Matkkials and Mktmods 

lei anther lobes growing downward into stiff, awn- 

, , 7 cv ! SAMPLING 

like horns; and (o) nectarv often a median knob in 


line with the floral axis. Cronquist (1981) referred Multiple taxa were chosen from each of the five 
to the Orobanchaceae as . . derived from Scro- major lineages ol parasitic: figworts and broomrapes 
phulariaceae,” but, due to its parasitic habit and found previously (dePamphilis et al., 1997). We 
parietal placentation, retained it as a separate lam- also included Schwalbea , a monotypic genus that 


ily. finally, a cladistic analysis of DNA sequences may represent a sixth parasitic lineage (Pennell, 
from the plastid gene rps2 showed the Orobancha- 1935). Both parasitic figwort tribes, Rhinantheae 
ceae to be firmly ensconced within the parasitic and Buchnereae, were well sampled. From the 
figworts, which are not monophyletic (dePamphilis 


broomrapes, we sampled Epifagus , ConopholLs , 
et al., 1997). Despite the lack of monophyly, we Boschniahia , Cistanche , and all four sections of 
use the name Serophulariaceae and its common Orohanche: Gymnocaulis ((). uniflora and 0. fasci- 

, figworts, as they have been traditionally culata ), Nothophyllon (O. corymbosa ), Orohanche 

(O. caryophyllaceae , 0. cernua , and 0 . hederae ), 
This study assessed the monophyly of three and Trionychon (0. ramosa). 

rps2 phylogeny showed that the nearest relatives 


used, until a reclassification is published. 


groups: genus Orohanche , the broomrapes, and tht j 


figw 


to the parasites are clearly not the Digitaleae 


the study also sought the sister group of each clade. (dePamphilis et al., 1997), as was previously 
The position of Schwalbea was also examined. thought. We therefore sampled members of all of 


Schwalbe a was thought to be among the most Wettstein’s (1897) figwort tribes except two small 
“primitive” of parasitic figworts, based on its fifth ones, Aptosimeae and Manuleae, for which material 
(posterior) sepal, septicidal dehiscence of the cap- was not available. Additional sampling outside the 
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/ k 


figworts was based 
(dePamphilis et al., 1997, and unpublished), and 
emphasized lineages that appeared to be closely 
related to the parasites. We also included Linden- 
bergia, which, although not parasitic, shares some 

rnilarities with the parasitic figworts (Briihl, 
1920). Trees were rooted with Ligustrum (Ole* 
and Nicotiana (Sola 
ed, along with collection information and the 
GenBank accession numbers of the sequences, in 

Table 1. 


) 

). The taxa used are list- 


DNA ISOLATION. AMPLIFICATION. AND SEQUENCING 

The CTAB method (Doyle & Doyle, 1987) was 
used to isolate total DNA from plant leaf or stem 


t 


. Amplification of rps'2 followed dePamphilis 
et al. (1997), and matK sequences were amplified 
using IX Taq Extender buffer (Stratagene), 0.2 mM 


dTTP 


da). 


3.6 mM MgCl 2 , 0.32 p,M each primer (see Fig. 1), 


0.25 unit of Taq DNA polymerase, 0.25 unit of Taq 


Extender coenzyme (Stratagene), and ~ 500 ng of 
total DNA in a 50 p.L volume. Figure 1 shows the 
primary matK primers (Genosys) used. Additional 
species-specific primers were occasionally needed, 
especially for Orobanche. PCR products were pu¬ 
rified in \% agarose gels using Qiaquick columns 
according to the manufacturer's instructions (Qia- 
gen). Sequences were generated by two methods: 
with the ABI 377 autosequencer (P. E. Biosystems) 

according to the manufacturer's instructions (with 


the exception that 


tions were done in 10 p.L, 


rather than the standard 20 |xL volume) and man¬ 


ually, with the Seq 


(U.S. Biochemicals) dou¬ 


ble-stranded method (dePamphilis et al., 1997). 
Both strands were sequenced, and all ambiguities 
clarified by individual examination using Se- 
quencher 3.0 software (GeneCodes, Ann Arbor, 
MI). Sequences were translated to verify that the 
protein-coding regions contained no internal stop 
codons, which would signify a possible sequencing 
error or pseudogene sequence. 


ALIGNMENT 

The rps‘2 alignment was simple, requiring only 
two small indels. The matK alignment was more 
complicated, requiring 14 locations with indels 
ranging from 3 to 21 bp, some of which varied 
among taxa in length and exact position. Initially 
22 alignments were produced, using Clustal W 1.4 
(Thompson et al., 1994) and the following varia¬ 
tions: gap opening penalties (GOP) ranging from 5 
to 30, gap extension penalties (GEP) from 1 to 10, 
with and without transition/transversion weighting. 


with the NJ guide tree replaced by a “conservative" 
parsimony tree (derived from a data set with length- 
variable regions removed), with complex length- 
variable regions coded as multistate characters, as 
in Baum et al. (1994), or with length variation di¬ 
vided into characters based on the longest possible 
independently varying units (Young, 1998), and, fi¬ 
nally, protein translations were aligned and the 
DNA alignments altered to match. 

The 22 alignments thus produced were evaluated 
according to the consistency among characters (the 
rescaled consistency (RC) index of Farris (1989)] 
on the resulting maximum parsimony (MP) trees 

(Wheeler, 1995). This provided an objective opti¬ 
mality criterion: whichever alignment yielded the 
MP tree(s) with the highest RC was considered the 
best alignment. 


PHYLOGENETIC ANALYSIS 


We used two methods of phylogenetic analysis in 
the program PAUP: maximum parsimony (MP) and 
maximum likelihood (ML). For rps2, positions ho¬ 
mologous to positions 48—660 of the Nicotiana gene 
were used [ Nicotiana plastic! genome positions 

16275—16887 (Shinozaki et al., 1986)], and for 
mat K, Nicotiana gene positions 297-1286 were 
used (Nicotiana plastic! genome positions 2425— 
3414). Because there were more than twice as many 
taxa sampled for rps2 than for matK, we analyzed 
each gene alone and in two different combined an¬ 
alyses: one with only those taxa sequenced for both 
genes, and one with all the taxa. In analyses of the 
raatK-only and small-merger matrices, heuristic 
searches were performed with 100 random addition 
orders and TBR branch swapping. The rps2-on\y 
and large-merger matrices were too large for a full 
heuristic search, so we used a search strategy de¬ 
signed to sample tree space thoroughly in large data 
sets (Catalan et al., 1997). In the mat K-oidy anal¬ 
ysis, the effect of the two regions with overlapping 
gaps was investigated by repeating the analysis 
with the gaps removed. Bootstrap support (Felsen- 
stein, 1985) was estimated from the mat K-only and 
small-merger matrices, using the same parameters 
as the heuristic search (Fig. 3), and at least 500 
replicates (Figs. 2, 4). For the r/w2-only and large- 
merger matrices, each bootstrap replicate was lim¬ 
ited to five random addition orders and five fully 
swapped trees each (Figs. 2, 4). Bremer support 
values (Bremer, 1988) were calculated for all four 
matrices. Topological constraints were used to find 
the number of extra steps that would be added to 
the MP trees under specific hypotheses (Table 2). 
Maximum-likelihood analyses were conducted on 
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the nm/K-only and small-merger matrices; the rps'2- 
only and large-merger matrices had too many taxa. 
The substitution model included a transition/trans¬ 
ratio of 2.0 and base frequencies estimated 
from the data (Hasegawa et al., 1985). Analyses 
were conducted using test version 4.0d54 of 
PAUP*, with permission of the programs author. 
David L. Swoflord. The aligned data can be ob¬ 
tained from TreeBASE (studv accession number 

w 

S402) at htti >://www.herl >aria. harvard, edu/treehase/ 


Rksilts 


Although a cpDNA phylogeny cannot he as- 
1 identical to the organisinal phylogeny (e.g., 
Doyle, 1992), introgression or lineage sorting are 
not likely to he problems above the level of genera, 
so we expect close correspondence between the 
cpDNA phylogenies and the organisinal phylogenv. 
The lour data matrices used in the phylogenetic 
analyses are detailed in Table 2. The first, referred 
to as 4 V/>.s2-only,” contains rps2 gene sequences 


from 63 taxa. The second matrix (mat K-only) con¬ 


tains mat K sequences from 26 of those taxa (plus 
Cistanche , just recently obtained). The third matrix 
(small-merger) combines sequences for the 26 taxa 
that have data from both genes (all of the taxa in 
the matK analysis except Cistanche , which had 


only two-thirds of its matK gene sequenced, due to 
technical difficult ies). This type of analysis (many 
characters, few taxa) is expected to provide the 
strongest support for clades on trees (Sanderson, 
1989). The fourth matrix (large-merger) also ex¬ 
cludes Cistanche , but contains all 63 remaining 
taxa; about 39% of the data are missing. This type 


of analysis (many taxa, fewer characters per taxon) 
can be expected to have lower support for clades 
on trees, but can give at least a preliminary indi¬ 
cation of relationships for all the taxa. 

The strict c 



s tree resulting from the rps2- 
onlv MP analysis (Fig. 2) supports the monophyly 
of the parasitic figworts and broomrapes and indi¬ 
cates that Lindenbergia is the nearest relative to the 


parasites, among those plants sampled. Within the 
parasites, early branching events are unresolved, 
but the transitional genera (Lathraea, Harveya, and 

w 

Hyobanehe) clearly do not group with the broom- 
rapes. Constraining them to do so adds 14 steps to 
the tree (Table 2). 

The best alignment of the matK sequences, re¬ 
sulting in MP trees with a RC — 0.380, was gen¬ 
erated with Clustal W alignment parameters GOP 
= 15, GEP = 5, transition/transversion weighting 

on, and the seven resulting indel characters were 


coded as in Baum et al. (1994). The indels, all in 
multiples of three bases, were distributed through¬ 
out the gene, but were more common near the ends. 
Because it had the highest RC, this alignment was 
chosen for the K-only, small-merger, and large- 
merger analyses. Similar alignment parameters pro¬ 
duced RC values only slightly lower and yielded 
the same phylogenetic* trees. However, parameters 
that differed substantially (such as GOP < 5 or 
GOP > 25) produced RC values less than 0.375 
and led to less resolved consensus trees. 

The mafK-only MP analysis resulted in six short¬ 
est trees, which differ only in the placement of 
Schwalbea and Cistanche . The strict consensus tree 
(fig. 3) resolves the basal branches of the parasitic 
clade better than the r/w2-only analysis. The matK- 
only analysis supports the monophyly of traditional 
Orobanchaceae, whose sister group is the Striga- 
Harveya clade (though support measures are not 
high). It also indicates that the genus Orobanche is 
diphyletic (polyphyletic, forming two clades), with 
moderately high support. It also resolves the Cas- 
tilleja clade as sister to the Melampyrum clade. The 
ML (maximum likelihood) tree differs only in that 
Antirrhinum and Veronica form a monophyletic 
group rather than a paraphyletic one. This indicates 
that the parsimony algorithm may have allowed Ve¬ 
ronica and Nicotiana to be attracted to each other 
due to their long branches (Felsenstein, 1978). That 
Antirrhinum belongs with Veronica is supported by 
both the rps2-on\y and large-merger analyses, 
which break up these branches with the addition 
of more taxa (Figs. 2, 4). When the two regions 
containing overlapping gaps were removed from the 
analysis, as well as the gap characters generated 
by these regions, the strict consensus MP tree dif¬ 
fered by a single feature: the sister relationship of 
two of the outgroups, Hemimeris and Verboseum, be¬ 


comes unresolved. 

The small-merger analysis (Fig. 3) gives addi¬ 
tional support to the findings of the previous two 


analyses. Cistanche has been removed, but other¬ 
wise the resulting MP trees have exactly the same 
topologies. Again, the ML analysis conflicts only in 
the placement of Veronica . The monophyly of the 
parasites (hemi- plus holo-) is strongly supported. 


as is the position ot Lindenbergia as sister group to 
the parasites. The holoparasites are clearly poly¬ 
phyletic. Although the broomrapes (as sampled) are 
monophyletic, Harvey a and Lathraea do not form a 
clade with them. Instead, each of these holoparas- 

a close relative of green hemi- 


ltes 


revealed 


parasites. Harveya is related to Alectra and Striga 
(and others of tribe Buchnereae; Fig. 4). Lathraea 
is in a clade with Tozzia and Melampyrum (and oth- 
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trnK 11 -53F 323F 775F 1210F 



510R 950R 1189R 1349R trnK 2621 


trnK 11 

trnK 2621 
matK -53F 
matK 323F 
matK 775F 
matK 1210F 
matK 1349R 
matK 1189R 
matK 950R 
matK 510R 


CTCAACGGTAGAGTACTCG 

AACTAGTCGGATGGAGTAG 

CTTGTTTTGRCTNTATCGCACTATG 

TATTCTCAAATCATATCRGAGGG 

TCTTGAACGAATCTATTTCTRYGG 

GGGCATCCCATTAGTAAGNC 

CTTTTGTGTTTCCGAGCYAAAGTTC 

CGGTTACTAATGGGATGCCC 

CCACARCGAAAAATRMCATTGCC 

GAAGAGTTTGAACCAAKAYTTCC 


tigure I. Map <>l the unit K gene (thick lint*) and Hanking regions (tliin line) with primer sites indicated. Primer 

sequences are listed 5' to .V: some contain lUPAC ambiguity codes. trnK I I was designed by G. beam; trnK 2921 
was designed by K. Steele (Steele & Vilgalys. 1994). 


ers of tribe Rhinantheae). Tozzia is therefore not an 


evolutionary transition to Lnthraea, though they are 
fairly close relatives. The genus Orobanche is di- 
phyletic (Fig. 3). Two N ew World species group 
with Boschniakia , and four Old World species 
group with Epifagus. The rbcL phylogeny also sup¬ 
ports Orobanche diphyly (Wolfe & dePamphilis, 
1998). However, the support values for the New 
World clade are low, so this result must he regarded 
as preliminary. 

The strict consensus tree resulting from the 
large-merger analysis (Fig. 4) is more resolved and 
generally better supported than that from the rps'2- 
only analysis, but is otherwise completely congru¬ 
ent with it. The relationships at poorly supported 
nodes, especially among the outgroups, should he 
regarded as very preliminary. Neither subfamily of 
Scrophulariaceae is monophyletic, and this is not 
only because the broomrapes, mints, and other fam¬ 


ilies derive from within the subfamilies. The trib 
Digitaleae s 



W^ettstein, which includes Veroni 
ca and Hemiphragma , forms a clade with the trib 
Antirrhineae and its relatives rather than with th 
rest of the subfamily Rhinanthoideae. Likewise nei 
ther parasitic tribe is monophyletic. Constraint an 
alyses show that many extra steps would have to b 
added to make the subfamilies and parasitic tribe 
monophyletic (Table 2). Surprisingly, there is stron, 
support for a clade that includes members of botl 
the Rhinantheae and the Buchnereae. All of th< 


genera in this clade have their center of distribution 
in North America (Mabberley, 1997). 


Discission 


These analyses of rps2 and matK gene sequences 
allow us to better understand the phylogeny and 
thus the evolution of parasitism within the figworts 
and broomrapes. The number of characters and 
taxa we have utilized exceeds that of previous stud¬ 
ies. Out of the 27 taxa contained in the rps2 tree 
and the matK tree, there are only two conflicts. The 
combination of the data sets is therefore supported 

(Mason-Gamer & Kellogg, 1996). The inonophyly 

of the parasitic figworts and broomrapes is strongly 
supported in all four analyses. The inonophyly of 
the parasites also indicates a single origin of par¬ 
asitism. It is not yet clear how many times haus- 
torial parasitism has evolved in other groups, but 
an estimate of at least 11 times (Nickrent et al., 
1998) suggests that the habit may be relatively un¬ 
common in flowering plants. However, once hemi- 
parasitism has been established, the loss of pho¬ 
tosynthesis (and evolution of holoparasitism) may 
be more common. 


Holoparasitism has arisen five independent times 
in the figwort/broomrape clade. In addition to the 
Harveya, Lathraea, and broomrape lineages, it has 
also occurred once within the genus Aleetra and 
once within the genus Striga (dePamphilis et al.. 
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Figure 2. Phylogeny of figworts and broom rapes based on rps2 gene sequences. The strict consensus tree is shown, 
with bootstrap values (above branch) based on 5(H) replicates and Bremer support values (below branch). Bold font 
typeface and bold font branch segments indicate holoparasitic taxa. 
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Alectra sessiliflora 

Harveya purpurea 

Striga asiatica 

Boschniakia hookeri 

Boschniakia strobilacea 

Orobanche corymbosa 
Orobanche uniflora 
Cistanche phelypaea 

Epifagus virginiana 
Orobanche caryophyllaceae 

Orobanche hederae 
Orobanche ramosa 

Orobanche cernua 

Castilleja liniariifolia 
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Lathraea clandestina 
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Kigelia pinnata 
Hemimeris sabulosa 
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Antirrhinum majus 

Veronica arvensis 

Nicotiana tabacum 


Figure* 3. Results of both the mat K-onlv analysis and the small-merger analysis. The topology shown is the strict 
consensus tree from tin* mat K-only analysis. Cistanche and Schwalbea have unresolved placements; when these taxa 
are removed, a single, fully resolved tree is found. The small-merger combined analysis of rps2 and null K yielded a 
strict consensus tree with this same topology, except that Cistanche was not included in that analysis. The four numbers 
displayed at each node represent, top to bottom, the mat K bootstrap value (based on 1500 replicates), the small-merger 
bootstrap value (based on 838 replicates), the mat K Bremer support value, and tin* small-merger Bremer support value. 
In one case, different species were used to represent a genus: Verbascum blattaria was sequenced lor rps2\ l. thapsis 
was sequenced for mat k. Bold lout typeface and bold ford branch segments indicate holoparasitie taxa. 
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f igure 4. Strict consensus tree from tlie large-merger analysis. Moots!rap values (above) based on 500 replicates 
and Bremer support values below. Bold font typeface and bold font branch segments indicate holoparasitic laxa. 
Subfamily codes: A = Antirrhinoideae. B = Rhinanlhoideae. Codes for Wettstein's tribes: U = Buehnereae. It = 


Bhinantheae. G 


Gratioleae. Ch 


Cheloneae, H 


Uemimeridae. V 


Verbasc 
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neae. I) = Digitaleae, Ca = Calceolarieae. Biogeography codes 
province: Am = Americas, Af = Africa, Ku = Eurasia. 


pneae, A = Antirrhi- 
taxa have their centers of distribution in the following 
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Table 2. Descriptive* measures from parsimony analyses and results of monophyly constraint 

analyses. Small = 

w 

small merger, Ixirge = 

large merger, MU = 

: maximum parsimony. 




rps2 

mat K 

Small 

I^arge 

INo. of taxa 

63 

27 

26 

63 

No. ol characters 

617 

12< >5 

1822 

1822 

No. of MU trees 

15000 

6 

3 

13941 

Free length 

727 

1421 

1823 

2137 

Cl 

0.552 

0.669 

0.677 

0.627 

RC 

Constraints 1 : 

0.370 

0.380 

0.398 

0.381 

K-S 

20 

21 

34 

36 

SF-A 

15 

15 

22 

29 

SF-R 

19 

23 

37 

38 

T-B 

6 

17 

22 

25 

T-R 

6 

17 

22 

22 

(;-() 

0 

9 

8 

7 

B&TG 

14 

54 

67 

67 

1 The number of extra steps added when 

the following groups are constrained to be monophyletk 

• (Swofford, l ( W3): 

K-S = family Serophui 

lariaeeae (excluding 

broomrapes). SF-A = subfamily Ant 

irrhinoideae, SF-R : 

= subfamily Rhin- 

anthoideae, T-B = tril 

>e Buchnereae, T-R 

= tribe Rhinantheae, G-0 = genus 

Orobanche , B&TG = 

= broomrapes with 

the transitional genera 

Harveya , hithraea , 

w 

and Hyobanche (if present on tree). 



1997). In all cases, holoparasitic taxa have been 


derived from hemiparasitic ancestors. This suggests 


that the loss of photosynthesis may only occur sub¬ 


sequent to the evolutionary origin of the haustori- 
um, the key character defining a parasitic strategy. 


THE EVOLUTIONARY SERIES IS NOT SUPPORTED 


The transitional genera Lathraea , Harvey a , and 
Hyobanche are each more closely related to green 
hemiparasitic lineages than they are to the broom- 
rapes. The placement of Lathraea with the Melam- 
pyrum group is found in 100% of bootstrap repli¬ 
cates and all trees up to 23 steps longer (Fig. 3). 
The grouping of Harveya (and presumably Hyob¬ 
anche, based on Fig. 4) with hemiparasites of the 
Striga group is also extremely well supported (boot¬ 
strap value = 100, Bremer support value = 19). 
This reveals that the characters used previously to 
group these taxa with the broomrapes (e.g., loss of 
green color, reduction ol leaves, the condensation 
ol the vegetative axis, and reduction ol cell number 
in the ovary) are all homoplastic and may be cor¬ 
related with holoparasitism. The past practice of 
grouping the transitional genera with the broom- 
rapes may be due to the fact that their specializa¬ 
tion involves the loss or reduction of many features, 
and convergently reduced or lacking features will 
seem to be homologies due to common ancestry. 


OTIIKH SYSTKMATIC CONCl.l SIONS 


Olmstead and Reeves (1995; Reeves & 01m- 
stead, 1998) began the disintegration of the Scro- 


phulariaceae by reporting that there are at least two 
major lineages of figworts. In studies of the chlo- 
roplast genes ndh F and rbcV, they identified the 
fact that Antirrhinum, Digitalis, and Veronica form 
a clade distinct from Scrophularia, Verbascum, and 
Selago. We have identified five additional clades of 
nonparasitic figworts. Leucophyllum represents a 
clade sister to the Myoporaceae. Mimulus and Pau- 
lownia represent clades that are not in the parasitic 
figwort clade, but are closer to it, Verbenaceae, and 
Bignoniaceae than they are to the Scrophularia 
clade. Finally, Calceolaria and Gratiola represent 
two clades whose positions are unresolved. Each of 
these lineages may end up as a family of its own 
or may be included with one of the other families, 
depending on its position on future, more resolved, 
phylogenies. A combined analysis is being con¬ 
ducted, using rps'2, ndhF, and rbcL and extensive 
sampling of nonparasitic Scrophulariaceae and re¬ 
lated lineages (R. G. Olmstead, C. W. dePamphilis, 
A. D. Wolfe, N. D. Young, W. J. Elisens & P. A. 
Reeves, unpublished). 

Based on Figure 4, we can see that neither of 
the traditional subfamilies (Antirrhinoideae, Rhin- 
anthoideae) are monophyletic. When trees are con¬ 
strained to contain a monophyletic subfamily, they 
are at least 15 steps longer (Table 2). In order to 
circumscribe monophyletic groups, the subfamilies 
must be broken up, and their members placed into 
at least five new or re-defined families, separated 
by other existing families (R. G. Olmstead, C. W. 
dePamphilis, A. D. Wolfe, N. D. Young, W. J. Eli- 
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sens & P. A. Reeves, unpublished). One redefined 
family will combine taxa out of both subfamilies. 

Digitalis and Veronica, classified by Wettstein 
(1897) in the Rhinanthoideae, are herein (Fig. 3) 
more closely related to Antirrhinum than to the re¬ 
mainder of Rhinanthoideae. This implies that the 
floral aestivation character used to distinguish 
Rhinanthoideae from the other subfamily, Antirrhi- 
noideae (Bentham, 1846; Bentham & Hooker, 
1876; Wettstein, 1897; Armstrong & Douglas, 
1989), has changed more than once in the evolution 
of the group: Digitalis, Veronica, and their relatives 
probably represent an independent origin of the 
rhinanthoid corolla aestivation type: the traditional 
Rhinanthoideae are clearly polyphyletic. This 
placement of Digitalis and Veronica agrees not only 
with ndhF/rbcL phylogeny (Olmstead & Reeves, 
1995) but also with pollen morphology. Pollen exine 
structure is tectate with reticulate sculpturing in 
Digitalis, Veronica, Antirrhinum, and Chelone, a 
structure that is otherwise rare in the family (Mink¬ 
in & Eshbaugh, 1989). 

The nearest relative to the parasitic figworts and 
broomrapes is the figwort genus Lindenbergia, not 
Gesneriaeeae. Gesneriaceae (here represented by 
Kohleria) must therefore have derived unilocular 
placentation independently from Orobanchaceae. 
Lindenbergia, a southeast Asian genus, has tradi¬ 
tionally been placed in the tribe Gratioleae in sub¬ 
family Antirrhinoideae (Bentham & Hooker, 1876) 
despite corolla aestivation similar to parasitic fig- 
worts and broomrapes (Cooke, 1903; Duthie, 1903— 
1920; Briihl, 1920; J. E. Armstrong, unpublished). 
Briihl (1920) has suggested that Lindenbergia 
shares a closer relationship with the parasitic fig- 
worts than with the Gratioleae, which agrees with 
our analysis (Figs. 2—4). Lindenbergia is here shown 
to be the sister group to the parasites, confirming 
its floral similarity. The clade containing all para¬ 
sitic figworts, broomrapes, and Lindenbergia is 
strongly supported. In the small-merger analysis, it 
has a bootstrap value of 99 and a Bremer support 
value of 11. In addition, this clade has a defining 
(though not unique) morphological synapomorphy: 
anterior lobes external in corolla aestivation. This 
group warrants family status, and the appropriate 
name is Orobanchaceae, a conserved name (Greu- 
ter et al., 1994). We recommend expanding the 
Orobanchaceae and suggest the following clade def¬ 
inition: 

(trobancliaceae are the least inclusive clade that con¬ 
tains Orobanc he uniflora, Schwalbea americana, and 

Lin den bergia philippine nsis. 

This definition, along with the changes suggested 
in Olmstead et al. (R. G. Olmstead, C. W. 


dePamphilis, A. D. Wolfe, N. D. Young, W. J. 
Elisens & P. A. Reeves, unpublished), forms a new 
taxonomic system, which is compared to Wettstein’s 
(Wettstein, 1897) in Table 3. Orobanchaceae now 
include hemiparasites and nonparasites, as well as 
holoparasites. 

Within the parasites, Schwalbea is among the 
basal lineages, as suggested by its morphology 
(Pennell, 1935), but its exact position remains un¬ 
resolved. The two parasitic tribes Buchnereae and 
Rhinantheae are each not monophyletic (Fig. 4). 
Agalinis (= Gerardia), Macranthera, and Seymeria, 
the New World representatives of Buchnereae 

Benth. 1846 (= Gerardieae Benth. & Hook., 1846), 

are part of a clade that is otherwise made up of 
genera of Rhinantheae with centers of distribution 
in North America. This clade is an example of the 
strong biogeographic pattern seen in this phylogeny 
(Fig. 4). Clades on the parasitic part of the phylog¬ 
eny tend to be made up of genera with their centers 
of distribution (Mabberley, 1997) in the same bio¬ 
geographic province (Fig. 4). There are large clades 
of taxa centered in Africa and Eurasia that include 
both hemiparasites and holoparasites. Even within 
the broomrapes a pattern emerges, with two clades 
of taxa centered in the Americas and one in Eur- 

QO|ll 

cir> i a • 

In Figures 3 and 4, the broomrapes are mono¬ 
phyletic. Their sister group (Fig. 3) is the Striga- 
Harveya group. However, this is not as well sup¬ 
ported as our other conclusions and must be 
regarded as preliminary. The other potential sister 
groups are the Tozzia-Lathraea and Castilleja-Sey- 
meria groups. Moderately well supported is the 
finding that Orobanche comprises at least two 
groups: one allied with Boschniahia, the other al¬ 
lied with Epifagus and Conopholis (Fig. 4). In the 
small-merger analysis (Fig. 3), these two new 
groupings are supported by bootstrap values of 70 
and 88, respectively. The current classification of 
Orobanche (Beck von Mannagetta, 1930; Collins, 
1973; Heckard & Chuang, 1975) contains two Old 
World sections and two New World sections. Beck 
von Mannagetta (1890) united the two New World 
sections into one branch of the genus, and the two 
Old World sections to form the other. Cytology in¬ 
dicates that one of the Old World sections, section 
Orobanche , has chromosome numbers that are gen¬ 
erally 2 n = 38, while the other three sections have 
chromosome numbers that are nearly always 2n = 

24, 48, 72, or 96 (Heckard & Chuang, 1975). To 

be consistent with Beck von Mannagettas scheme, 
the cytology suggests that section Orobanche might 
be monophyletic and derived from the other Old 
World section, Trionychon (including O. ranwsa, 2 n 
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Table 3. Genera of the Seroplmlariaceae/Orobanchaceae complex used in this study, deposited according to the 
taxonomic systems of Wcttstein (I8 ( >7), and Olmstead et al. (R. G. Olmstead, C. W. dePamphilis, A. 1). Wolfe, N. 1). 
Young, W. J. Klisens & P. A. Reeves, unpublished). 

Wettstein 

Olmstead et al. 


Schrophulariaeeae 

Scrophulariaceae 


Pseudosolaneae 

Verbose um 


Verbasceae 

hmcophyllum 


Verbose um 

Scrophularia 


Leucophyllum 

Selago 


Antirrhinoideae 

Hemimeridae 

Antirrhinaceae 


Hemimeris 

Antirrhinum 


Calceolarieae 

Chelone 


Calceolaria 

Gratiola 


Antirrhineae 

Hemiphragma 


Antirrhinum 

Veronica 


(dieloneae 

Digitalis 


Scrophularia 

Cal lit riche 


Chelone 

Hippuris 


Paulownia 

(iratioleae 

Calceolariaceae 


Mini ulus 

Calceolaria 


Lindenbergia 

Gratiola 

Orobanchaceae 


Selagineae 

Lindenbergia 


Selago 

Melasma 


Rhiuunthoidene 

Alert ra 


Digitaleae 

Macranthera 


Hemiphragma 

Seymeria 


Veron ica 

Agalinis 


Digitalis 

Sopubia 


Gerard ieae 

Huchnera 


Melasma 

Cycnium 

w 


Alert ra 

Striga 


Macranthera 

Haney a 


Seymeria 

w 

Hyobanche 

w 


A gal inis 

Castilleja 


Sopubia 

()rthocarpus 


Huchnera 

Triphysaria 


Cycnium 

✓ 

Melampyrum 


Striga 

Tozzia 


Han'eya 

hithraea 


Hyobanche 

m 

Euphrasia 


Rhinantheae 

Parentucellia 


Castilleja 

Barts ia 


()rthocarpus 

Pedicular is 


Triphysaria 

Hhinanthus 


Melampyrum 

Iximourouxia 


Tozzia 

Schwalbea 


Ixithraea 

Cistanche 


Euphrasia 

Conopholis 


Parent ucellia 

Epifagus 


Hart si a 

Hoschniakia 
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Table 3. Continued. 

Wettstein 

Olmstead et al. 

Pedicular is 

Orobanche 

Rli inant has 


Ijamourouxia 

No family designation, near Orobanchaceae 

Sell walbea 

Mini ulus 

Orobanehacae 

Paulownia 

Cistanche 


Conopholis 

Bignoniaceae 

Epifagus 

Kigelia 

Boschniakia 


Orobanche 

No family designation, near Bignoniaceae 

Bignoniaceae 

Schlegelia 

Kigelia 


Sc hiegelia 

Verbenaceae 

Callitrichaceae 

Verbena 

Callitriche 


Gesneriaccae 

Gesneriaceae 

Kohler ia 

Kohleria 

Hippuriaeeae 


Hippuris 


Myoporiaceae 


Myoporum 


Verbenaceae 


Verbena 



= 24). Our plastid phylogeny (Fig. 4) supports the 
Old World—New World division, hut finds a para- 
phyletic relationship of section Orobanche ((). car- 
yophyllaceae, 0. hederae, and 0. cernua , all '2n — 
38) to section Trionychon (2n = 24). This indicates 
that the Old World branch of the genus may have 
begun with 2 n = 38 and later experienced a re- 


ducti 


ion in 


k'V' 


tion Trionychon to 2 n 


24. No 


counts have been reported lor Epifagus , Conophol- 
is , or Boschniakia. Cistanche has 2 n — 40 (Ham- 

bler, 1956). 

Having a classification that reflects monophyletic 
relationships will be a great advantage, especially 
to comparative biologists. These phylogenetically 


defined groups may even be 


to identify based 


on morphological characters, compared to the tra¬ 
ditional family Scrophulariaceae, which is recog¬ 
nized by symplesiomorphies such as the presence 
of endosperm, capsular fruit, and strongly zygo- 
morphic flowers. Furthermore, Antirrhinum majus 
serves as a model organism in the field of devel¬ 
opmental biology (Coen & Nugent, 1994; Bradley 
et al., 1996). Such work can be placed in an illu¬ 
minating context by the study of its relatives and 


their phylogeny (R 


& Olmstead, 1998). 


Knowing that photosynthesis has been lost mul¬ 
tiple times in Orobanchaceae opens up many op¬ 
portunities for comparative analysis. It allows for 
the comparison of rates of DNA base substitution 


(dePamphilis et al., 1997), plastid genome struc¬ 
tural evolution (dePamphilis, 1995), as well as rates 
of loss of particular photosynthetic genes (Wolfe & 
dePamphilis, 1997). Haustorial anatomy and phys¬ 
iology, host plant use, morphology, as well as the 
genetic changes that have accompanied parasitism, 
can all be investigated using the phylogeny for ref¬ 
erence. Awareness of the multiple, independent or¬ 
igins of holoparasitism provides us with a powerful 
comparative framework in which to study the pro¬ 
cess of parasitic evolution. 
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